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1 A Recursion for the Farey Sequence 

Represent a rational number / as an ordered pair (n, d) where n is the numerator 
and d is the denominator. The Farey sequence of order m, F^, is given by the 



recursion 

' 'i-i + m 



i+l 



/. - /, 



j-1 



di 

where /i = (0,1) and /a = (l,m) f |Har02j . |Farl6j . IHalTOQ . 

The Farey sequence Fm is the sequence of irreducible rational numbers {^} 
with < n < d < m arranged in increasing order. It is occasionally useful to 
have available an expression for -Fm+i in terms of Fm- 

2 A Recursion for the Farey Sequence Sequence 

Represent an element / of the Farey sequence F„i as an ordered triple {nf,df,s) 
where n/ is the numerator, is the denominator and s is the number of Farey 
sequences following Fm until a new fraction is created immediately after /. 
When the context is clear, the subscript / will be omitted. 
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In this notation F2 is given by 



^^2 = (0,1,1)11(1,2,1)11(1,1,0). 

Starting with the null sequence, Fm+i may be computed from as follows: 
For i from 1 to \ F^ \ —1 

If{Si>l) 

Fm+l <= Fm+1 \\{ni, di, Si - 1) 

else 

^ Fjn+i\\{ni, di, di)\\{ni + ni+i,di + dj+i, dj+i) 

Next i 

F„+i ^F„+i||(l,l,0) 

We shall refer to the second (mediant) fraction in the else clause as a created 
fraction and denote the set of all created fractions in the transition from F^ to 
Fm+i by Cm- The number of elements in Cm is <^(m + 1) where is the Euler 
totient function. We shall denote the (initial) value of s of a created fraction / 
by Sf. 

Note that if (n^, d, 1) G Fm for 1 < d < to then \Cm\ = firn + 1) = m and 
m + 1 is a prime. 

3 Properties 

The following are some elementary properties of the Farey sequences generated 
using the above notation: 

(1) Fm is the Farey sequence of order m 

^ (n, d, s) e Cm 'i^d = m+ l,n<d and n\ d. 

(2) {rii, d, Si), {nj,d, sj) e Fm and j ^ Sii^ sj. 



(3) (n, d,Sf) e Cm' and (n, d, s) e Fm for m > m' + s/ 

^ (n, d, s) e F„+fed for A; = 1, 2, . . . 

(4) / = (n, d, s) G =^ s = d — (m — s/) mod d 

(5) / = (n, rf, 1) G Cm <^ (m - s/ + 1) mod d = 

(6) (rii, di, 1), {nj,dj, 1) G with di^ dj 

(7) m + 1 is prime (rij, i, 1) G for i = 1, . . . , m 
These properties may be interpreted as follows: 

(1) says that if Fm is the Farey sequence of order m then F^+i is the Farey 
sequence of order m + 1. 

(2) says that all fractions with the same denominator in a given Farey se- 
quence have different s values and in particular the s / values of elements of Cm 
all different. 

(3) says that after an initial period of length Sf the s value associated with 
a fraction with denominator d cycles with a period of length d. 

(4) is a restatement of (3). 

(5) is a special case of (4). 

(6) says that the denominators of the fractions immediately preceding cre- 
ated fractions are all different. As the sum of the preceding and succeeding 
denominators of a created fraction is m + 1, the same holds for the succeeding 
denominators. 

(7) is the special case of (6) where w + 1 is a prime. 

4 Functions of the Farey Sequence Sequence 

The above recursion for generating the sequence of Farey sequences may be 
used to describe properties of Farey sequences as a function of the order of the 
sequence. 
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For example, the distance between a Farey fraction / = (n, d, s) in Fm and 
the next larger Farey fraction in Fm as a function of m is given by 



As another example, one can express the order index of an element / of the 
Farey sequence Fm as a function of Fm and m as follows. Let i / denote the 
(initial) order index of / in Fd^. Then the order index of / in Fm for m > df 
is given by 



'm-Sg' 




'df-Sg 






. dg _ 



The result of Franel and Landau |FL24j regarding the Ricmann Hypothesis 
and the Farey series may be written using the above notation as 



The above recursion for generating the sequence of Farey sequences may also 
be used to generate sequences of primes. For this purpose, we develop notation 
to describe the cycling behavior of s. 





where $ is the Euler summatory function. 



5 Cycles 



If 



"^c(/, k) = < 



kdf + Sf — c 



C<Sf 



{k + l)df + Sf — c c > Sf 



then 



M,{f)^{m,{f,k),k = 1,2,...} 
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is the set of all m such that the value of s associated / in Fm is c. Further, if 

Mc{d) = U M,{f), 

df=d 

then M.c{d) is the set of all m such that the value of s associated with a Farey 
fraction in Fm with denominator d is c. The M.c{d) sets can be to define prime 
sieves and prime recursions. 



6 Prime Sieves 

p is a prime if 

p-i 

p-ie {^MM 
d=i 

This computation is equivalent to ensuring that there exists an f ^ Cd such 
that 

(^)- 

for each 1 < d < p — 1. 

p is the lesser of a twin prime |Guy04| if 

p-i 

P-lef] [Miid)f]M3id)} 

d=l 

p is a prime since 

p-i 

p-le f]Miid). 

d=l 

Further since {[p/2],p,2) G Fp we have {[p/2],p,l) e Fp+i along with (0,1,1) 
and + 1, 1). The condition 

p-i 

p-le f]M3{d) 
d=i 
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ensures that (n^, d, 1) e Jp+i for 2 < < p — 1. Therefore (n^, d, 1) e -Fp+i for 
1 < < p + 1 and hence p + 2 is a prime. 

7 Prime Recursions 

Let ma = Mi{2) and ds = 3. Set 

mj+i = mjPl Ali(dj) 

and 

di+i = 1 + niinm,+i. 

Then is the i*^ prime. 

Let (p, p + 2) be a twin prime and let 



q = l + mmf]J^Mi{d)f]Ms{d)] 



Then {q, q + 2) is the next twin prime. This recursion can be expressed incre- 
mentally as follows. 

Set mi = Mi{3)f]M3{3) and di = 3. Let 

di+i = 1 + min mi 



and 

di+i 



mi+i=mif]\ fl Mi{d)f]M3id) 

d=d, + l 



Then {di, di + 2} is the « prime pair. 
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8 Twin Prime Generation 



The following is a Mathematica module to generate the sequence of twin primes. 

TwinPrimes [twins_, k_] := Module [{i, j, 1, p, q}, 
1 = Intersection [Ems [3, 1, k] , Ems [3, 3, k]]; 
P = 3; 

q = Min[l] + 1; 

For[i = 1, i < twins - 1, i++, 

Print ["Twin Pair #", i, ": ", {p, p + 2}] ; 
For[j = p + 1, j <= q, j++, 
1 = Intersection [1, Ems[j,l,k], Ems[j,3,k]]; 

]; 

p = q; 

q = Min[l] + 1; 
] 

Print["Twin Pair #", i, ": ", {p, p + 2}] ; 
Print ["Twin Pair #" , i + 1 , " : " , {q, q + 2>] ; 
] 

where Ems [d , c , k] is 

y {m,if,k),k = l,2,...,k} 

df=d 

9 Software 

A Mathematica notebook implementing all the above computations is available 
from the author upon request. 
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